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Abstract 



The vortex system in a high- Tc superconductor has been studied numer- 
ically using the mapping to 2D bosons and the path-integral Monte Carlo 
method. We find a single first-order transition from an Abrikosov lattice to 
an entangled vortex liquid. The transition is characterized by an entropy 
jump /S.S ~ 0.4 kg per vortex and layer (parameters for YBCO) and a Lin- 
demann number cl ~ 0.25. The increase in density at melting is given by 
A/9 = 6.0 X 10~^/ A(T)^ . The vortex liquid corresponds to a bosonic superfluid, 
with ps = p even in the limit A — > oo. 
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Our understanding of the phase diagram of type II superconductors has improved sig- 
nificantly since the mixed state was introduced by Abrikosov in 1957 The vortex state 
is of particular interest for the high-T^ superconductors, where strong thermal fluctuations 
lead to melting of the vortex lattice and the appearance of a vortex liquid phase. Exper- 
imental evidence for a first order vortex lattice melting transition has been obtained from 
the observation of jumps in the resistivity and in the magnetization More recently, 

the latent heat of the transition has been measured directly in an untwinned YBCO single 
crystal 0. Theoretically, vortex lattice melting has been studied using various approximate 
techniques, including the renormalization group, perturbative expansions, density functional 
theory, and the Lindemann criterion 0. The absence of one simple and reliable theory has 
provoked large interest in numerical simulations. A number of models, such as the 3D frus- 
trated XY-model, the lattice London model and the lowest Landau level approximation, 
have been used, with no consistent picture emerging, however [0. A careful analysis of the 
3D vortex system has been carried out by Sasik and Stroud using the Lowest Landau 
Level approximation at constant applied field. They obtained a jump in the magnetization 
and could accurately trace the first-order melting transition in YBCO. 

A very fruitful concept was introduced by Nelson 0, who showed that the classical 
statistical mechanics of the vortex system can be mapped onto the quantum statistical 
mechanics of a 2D system of Yukawa bosons, i. e., bosons interacting with the potential 
V{R) = (7^Ko(-R/A), where Kq is a modified Bessel function, A is the London penetration 
depth, and is a coupling constant to be defined below. We will refer to this system as 
the Bose model of the vortex system. A consequence of this mapping is the prediction of a 
melting transition into an entangled vortex liquid, corresponding to a lattice to superfiuid 
transition in the Bose system. Unfortunately, perturbation theory does not work for strongly 
interacting 2D bosons and the only quantitative results are from ground-state Monte Carlo 



simulations []I0[, which have not studied the superfiuid density or the excitation spectrum. 
In this letter, we present the first path-integral Monte Carlo results for 2D Yukawa bosons. 
Both the cases A finite and A = oo (2D Bose Coulomb Liquid) are considered. The results 
are interpreted for the solid-liquid transition in the vortex system. 

In the Feynman path integral formulation, a d-dimensional system of massive quantum 
particles is equivalent to a classical d + 1-dimensional system of interacting elastic strings. 
The dimensionless imaginary-time action for N Yukawa bosons is given by 



'5A=M- Ettt^FE +EKoR^ , (1) 




where all energies are measured in units of g^ and all lengths in units of the particle distance 
Oo in the crystal, Og = 2/p\/3, with p being the density. The de Boer parameter A, A^ = 
/ma^g"^, measures the size of quantum fluctuations, (3 = g^ jT^ is the inverse temperature, 
and the Rj denote the particle positions. The partition function is the sum over all world 
lines weighted by this action and subject to the boundary conditions Ri(/5) = Rj(0), i. e., 
every line ends on itself or on some other line. This action is also valid in the limit A — > cxo 
if a uniform background charge is subtracted. In this case, the Bessel function reduces to a 
logarithm. 

Redefining parameters, Eq. (|l|) can be interpreted as the free energy for the vortex system 
in a type II superconductor, S /Ti = T jT [||, with 
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A = ^^, /? = ^, (2) 
ao^/2eleo T 

where T is the temperature of the vortex system and is the thickness of the sample. 
For an anisotropic superconductor, the elasticity is ei ~ s'^Eq, where = m/M < 1 is the 
anisotropy parameter. Two approximations are required for this mapping from bosons to 
vortices: First, the original London functional for the free energy contains retarded and 
advanced interactions between the vortex lines, which are mediated by gauge fields in the 



Bose picture [^. Here, retardation is neglected, equivalent to keeping only the first term 



in an expansion around straight lines. Second, it has been argued that vortex loops in the 



a6-planes are important for vortex lattice melting [|1^ . A simple estimate for the free energy 
of a loop of length L, T ^ {L/^){eeQ^ — Tin 3), shows that loops proliferate in the critical 
regime close to Hc2- We do not consider the critical regime in this work. 

The phase diagram for the Bose model, shown schematically in Fig. [I|, contains three 
phases: A classical high temperature normal liquid phase, a crystal for low temperatures 
and small quantum effects, and a superfiuid as quantum effects start to dominate at low 
temperatures. This can be understood by considering the three energy scales involved: The 
transition from a normal liquid to a lattice is determined by the competition between the 
thermal energy = g"^ / (3 and the potential energy g^. In the limit A — oo and A = 0, the 
transition takes place at j3m ~ 140 . With increasing quantum effects, we find a transition 



from a normal liquid to a superfiuid when the thermal energy matches the kinetic energy, 
A^/3 K, 1. At low temperatures, the competition between potential and kinetic energies 
determines whether the system is a crystal or a superfiuid. For /?, A = oo, it is known that 

A^ = 0.062 [|ig]. 



The boson phase diagram can be reinterpreted in terms of the vortex system, where (3 is 
proportional to the sample thickness and A measures the strength of thermal fluctuations. 
For thin samples, (5 < /9m;l/A^, we find a disentangled vortex liquid. In thicker samples, 
the system is either a lattice or an entangled vortex liquid depending on temperature and 
magnetic field. Note the non-trivial mapping between the H — T and A — /3 phase diagrams: 
In a thin sample, the constant field line (dash-dotted in Fig. |^) passes through the crystalline 
phase (low T), the disentangled liquid phase (intermediate T), and the entangled liquid phase 
(high T). With increasing L^, this line moves to higher values of /3 and the vortex lattice 
melting line is determined solely by the value of A^. 

The simulations are carried out using the path integral Monte Carlo technique, which 
is exact for bosons ||T^. Typical runs involve = 64 world lines, though systems as large 
as A^ = 100 have been used for analyzing finite-size effects. With (3 = 300 we accuarately 
capture the ground-state behavior, e. g., we reproduce the T-^ = result of Ref. [|1^] to within 
less than 1%. We find that M = 100 Trotter slices are sufficient to eliminate systematic 
errors in the bosonic quantum phase transition, see Ref. []T^ for details. For the largest 



systems we used roughly 30000 sweeps to equilibrate and 80000 sweeps to measure. 

The lattice to liquid transition is identified through the vanishing of the first Bragg peak 
(Q = Qi, r = 0) in the structure factor 

5(Q,r) = l(pQ(rKQ(0)), (3) 
where Pq(t) is the partial Fourier transform of the density operator, p(R, r) = 
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X]j<5 [R — R.i(T)]- The superfluid density of the Bose system is measured using the winding 
number [[14| , 

ei^^ (4) 

with the winding vector W = J2i lo '^'^ "^^R-j measuring the diffusion of the center-of-mass 
of the system in imaginary time. As the equihbration of the winding number is very slow 
for large systems, we use systems with = 36 to compute the superfluid density. For larger 
systems, we define the parameter pe according to 

Pe entangled lines 



p total number of lines 



where a line is entangled if it does not end on itself, Ri(/5) 7^ Rj(0). The parameter pe 
measures the importance of quantum effects in the system, and is therefore related to super- 
fluidity. We emphasize, however, that pe is not the superfluid density. 

We begin by considering the incompressible case with \ = 00. In Fig. 2 we show the 
results for the first Bragg peak and the parameter pe- The lattice disappears in a sharp 
transition at A^ ~ 0.062, in perfect agreement with ground-state simulations of the same 
model |]IU| . The height of the Bragg peak is related to the Lindemann number according to 
S{Qi) = A^exp(— Svr^cl/S), and we find cl ~ 0.25. Using Eq. (|^), we obtain the melting 
line 

B^(T) = 1^^. (6) 



as expected from a Lindemann criterion with A^ oc c| |TT|], giving good agreement with 
experimental results on YBCO. Its universal character applies only to large magnetic fields; 
the finiteness of A leads to a reentrant melting line at low fields . 

If the vortex lattice melts in a single transition, symmetry requires it to be first order 
16| . In Fig. we plot the energy per line and unit length e = {J-')/NLz. The corresponding 



energy of the vortex system is defined as = T'^Ot InZ/NLz, with Z the partition function 
|T7| . Using the scaling form JF = eoaof[{Rij/ao}], valid for large A, we obtain = e(l + 
t^) / (1— t^), where t = T/T^. The jump in entropy at the transition is given by T^^s^ = Ae^. 
From Fig. 3 we obtain Ae ~ O.OlSeo and therefore 

As^fcs /length] ^ 0.03£o(0)/T„. (7) 

Using parameters for YBCO (layer separation d = 12 A, Xab ~ 1400 A, and ~ T^) we 
obtain As^ ~ OAkg per vortex and layer, which compares favorably with the experimental 
result of As0 ^ 0.45 fc^ §. 

For the compressible case with A finite, the statistical attraction between the bosons 
produces an increase of the density upon melting. This maps to a densified vortex liquid 
due to the entanglement of the flux lines |T^. In order to study this effect, we have developed 
an isobaric Monte Carlo algorithm which allows us to fix the external pressure and let the 
volume adjust |1^. The results are shown in Fig. ^ for a system with A ~ 1.06ao. The 
transition is shifted to a slightly smaller value of A^ due to the weaker interaction between 
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the lines. The small shift, less than 4%, shows that Eq. (Q) gives a good description of 
the melting for a large range of fields. The change in density depends on the value of 
A/oq; in the present case Ap/p ^ 0.0003, which is consistent with the result (0) via the 
Clausius-Clapeyron relation, Ae/eo = 87rApA(T)^. 

We turn to the discussion of the liquid phase. In Fig. |^ we show the entanglement 
parameter pe which rises sharply at the transition, indicating that the vortex liquid entan- 
gles immediately upon melting. This result is in agreement with recent flux transformer 
experiments, showing that the vortex correlation along the magnetic field disappears at the 
melting transition [2U|. The possibility of a disentangled vortex liquid has attracted much 



interest in recent years ||2T[]. Two theoretical arguments have been put forward in favor of 
the existence of this phase: To begin with, the melting of the vortex lattice into an entangled 
vortex liquid involves the change of two symmetries: The transverse translational symmetry 
of the lattice and the longitudinal gauge symmetry. If these two symmetries do not change 
simultaneously, an intermediate phase will appear. A second argument stems from the anal- 
ysis of the 2D Bose Coulomb liquid: It can be shown that the suppression of long wavelength 
density fluctuations in this system leads to enhanced phase fluctuations and an algebraic 
decay of the off-diagonal long-range order (ODLRO) even in the ground-state (T^ = 0) 
[ pHI . However, the absence of a = Bose condensate has no straightforward implication 
for the superfluid density, which is related to the excitation spectrum rather than to the 
ground-state properties. In the inset of Fig. 2 we show the superfluid density measured by 
the winding number. Apart from a slight broadening of the transition due to the smaller 
system size, the result shows that ps = p as soon as the translational symmetry is restored 
in the liquid. Retardation may modify this result in the following ways: i) The decrease in 
the effective mass of the bosons (elastic tension of the vortices) favors the entangled state, 
ii) The retardation may render the entanglement unstable, thereby favoring a disentangled 
liquid. Thus, the question regarding the possibility of a disentangled liquid phase in the 
retarded model has not been completely settled. 

Additional information on the properties of the Bose superfluid/vortex liquid is provided 
through the analysis of the dynamic structure factor. Following Nelson the partial 
Fourier transform S{Q,t) takes the form S{Q,t) ^ S{Q,0) exp{— e(Q)|r|}, where e{Q) is 
the excitation spectrum of the Bose system. Thus, the bosonic excitation spectrum deflnes 
a longitudinal correlation length Ir = T/2A,.eo in the vortex fluid, where A^ denotes the 
roton minimum. We compute the excitation spectrum from our simulations by fltting the 
measured S{Q, r) to the single mode approximation, 

where e{Q) and T{Q) are the energy and inverse lifetime of the excitations. In Fig. ^ we show 
the resulting spectra both for the incompressible (A = oo) and the compressible (A < oo) 
fluids. Most interestingly, the phonon branch at small Q turns into a plasmon branch as 
A — >■ oo, while the roton minimum undergoes no visible change. From the roton minimum 
Ar ~ 0.027 we determine the correlation or entanglement length at the melting transition. 



Ir ~ l.GaQy Ei/eo, independent of the interaction range A. Note that in our simulations the 
roton gap is much larger than the temperature and we expect to probe the ground state 
behavior of the system. 



5 



In conclusion, our simulation of the 2D Coulomb Bose model reveals a single ^ 
quantum phase transition from a crystal to a superfiuid phase at ~ 0.062. This translates 
to a first-order melting transition of the Abrikosov vortex lattice into an entangled vortex 
liquid phase. The long-range interaction changes the bosonic excitation spectrum from 
phonons to plasmons for small Q, but does not modify the roton minimum. We find that 
superfiuidity is stable against static long-range interactions. 

We are grateful to M. Dodgson, M. Feigel'man, V. Geshkenbein, L. loffe, A. van Otterlo, 
and H. Tsunetsugu for interesting discussions. One of us (HN) is particularly grateful to 
H. Tsunetsugu for ideas on how to treat the long-range interaction, and to D. Ceperley for 
valuable advice on the algorithm. We acknowledge financial support from the Swiss National 
Science Foundation. 
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FIGURES 
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FIG. 1. Schematic phase diagram for a system of 2D charged bosons in terms of A and p. The 
solid lines represent phase transitions and quantum effects are relevant in the shaded region. In 
the vortex system, the parameters map to = T^/2e;eoflo ^^'^ ~ "^^o^z/T. The constant field 
line for a thin sample {Lz < JOT/eq) is shown (dash-dotted) as it runs through all three phases. 
For thicker samples, this line is pushed upwards. 
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FIG. 2. The first Bragg peak and the entanglement parameter for a system with 64 lines 
and (5 = 300. A sharp transition from a crystal to an entangled liquid is found at Am = 0.062. 
The structure factors just before and after the melting transition are displayed. The inset shows 
the superfluid density for a system with 36 lines. 
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FIG. 3. The energy e per line of the vortex system (right axis) for A = oo (sohd hne) and 
A 1.06ao (dotted hne). The energy of a perfect lattice with the same density has been subtracted 
in both cases. In the compressible system the transition shows a jump in the density (dashed line). 
Note the (small) shift in as the interaction range is reduced. 




FIG. 4. The excitation spectra £{Q) in units of g'^ for a system of 64 lines. Increasing the range 
A of the interaction shifts the sound mode to the plasma frequency but leaves the roton minimum 
unchanged. The latter collapses upon crystallization. 
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